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ABSTRACT 

We present a measurement of faint source confusion in deep, wide- field 4 /im 
images. The 1.8° x 1.8° images with 17" resolution are centered about the nearby 
edge-on spiral galaxies NGC 4565 and NGC 5907. After removing statistical noise 
and gain fluctuations in the focal plane array, we measure spatial fluctuations 
in the sky brightness to be bvl v = 2.74 ± 0.14 nW mT 2 sr _1 , approximately 
1% of the diffuse background level observed in a single pixel. The brightness 
fluctuations are confirmed to be associated with the sky by subtracting sequential 
images of the same region. An auto-correlation analysis shows the fluctuations 
are well described by unresolved point sources. We see no evidence for surface 
brightness fluctuations on larger angular scales ( 2' < 9 < 28' ). The statistical 
distribution of brightness fluctuations allows us to estimate the density of sources 
below the detection limit. We present a Monte-Carlo analysis of the undetected 
sources, which appear to be dominated by galaxies. Galaxies producing detected 
sky fluctuations contribute vl v {> S) = 1. 04^34 nW vrT 2 sr _1 to the cosmic 
infrared background, evaluated at S = 4.0 x 10~ 8 nW wT 2 . From the fluctuation 
data we can determine the integrated source counts N(> S) = I^Iq^q x 10 7 
sr _1 , evaluated at S — 4.0 x 10~ 8 nW m~ 2 . The observed fluctuations are 
consistent with reddened i-Cband galaxy number counts. The number counts of 
extracted point sources with flux vF u > 6.3 x 10~ 7 nW mT 2 are dominated by 
stars and agree well with the Galactic stellar model of Wright & Reese (2000). 
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Removing the stellar contribution from DIRBE maps with zodiacal subtraction 
results in a residual brightness of 14.0 ± 2.6 (22.2 ± 5.9) nW m~ 2 sr" 1 at 3.5 
(4.9) jum for the NGC 5907 field, and 24.0 ± 2.7 (36.8 ± 6.0) nW m~ 2 sr' 1 at 
3.5 (4.9) /im for the NGC 4565 field. The NGC 5907 residuals are consistent 
with tentative detections of the infrared background reported by Dwek & Arendt 
(1998), Wright & Reese (2000), and Gorjian, Wright & Chary (2000). 

Subject headings: cosmology: observations - - diffuse radiation - - galaxies: 
individual (NGC 4565, NGC 5907) — galaxies: statistical — infrared: galaxies 

1. Introduction 

Measurements of the cosmic infrared background (CIRB) probe the integrated emission 
from distant galaxies (Hauser and Dwek 2001). The brightness and spectrum of the CIRB is 
sensitive to galaxy number density, luminosity and evolution, as well as cosmology (Partridge 
and Peebles 1967; Kaufman 1976; Stabell and Wesson 1980; Dwek and Arendt 1998; Gispert 
et al. 2000). However, a search for the CIRB must account for the local zodiacal and stellar 
foregrounds. These have proven difficult to eliminate, even at A = 3 — 5 /xm where the local 
foreground from thermal emission and scattering by zodiacal dust is minimal. 

DIRBE reported upper limits on the CIRB at near-infrared wavelengths (Hauser et 
al. 1998), limited by foreground subtraction. The total astrophysical emission observed by 
DIRBE in a 0.7° beam contains a significant emission component from stars. Dwek & Arendt 
(1998) used ground-based i^-band galaxy counts along with the DIRBE data to obtain a 
tentative detection of the CIRB. Gorjian, Wright & Chary (2000) independently claimed a 
tentative detection of the CIRB using ground-based L-band counts to account for stellar 
emission in the DIRBE beams. Wright & Reese (2000) reported a detection using a different 
zodiacal model than used by Hauser et al. (1998). Matsumoto (2000) analyzed observations 
with the Near InfraRed Spectrometer (NIRS) on IRTS, with an 8' x 8' field of view that is 
significantly smaller than that of DIRBE, and reported an isotropic near-infrared component 
peaking at ~ 1.5 fim. 

Brightness fluctuations also probe the CIRB. The fluctuations contain information on 
the number counts and clustering properties of sources below the detection limit (Bond, 
Carr, and Hogan 1986; Fabbri 1988; Cole, Treyer, and Silk 1992; Kashlinsky and Odenwald 
2000). Poissonian fluctuations ("confusion noise") are predicted based on models of galaxy 
evolution, including stellar and cirrus components from our galaxy (e.g. Franceschini et al. 
1991). The brightness distribution depends on the index, a, of integrated source counts N(> 
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S) = N (S/So)~ a , and can thus be used to probe the density of sources below the detection 
limit (Condon 1974). The spatial distribution of fluctuations measures the clustering of 
sources (Knox et al. 2001) on size scales of galaxy clustering. For example, NIRS/IRTS 
observations show brightness fluctuations prominent at A6 m 100', appearing at the level 
of ~ 5% of the total diffuse sky brightness (Matsumoto 2000). Confusion noise sets the 
limiting sensitivity for future space-borne infrared observations, such as for SIRTF (Fanson 
et al. 1998) and Astro-F (Onaka et al. 2000). 

We present integrated source counts and a statistical analysis of residual brightness 
fluctuations in wide-field images obtained with a rocket-borne near-infrared camera at A = 4 
/im. These images provide a direct account of the point source contribution to the near- 
infrared sky brightness, and an indirect measure of the brightness contribution from sources 
below this detection limit. 



2. Observations 

The near-infrared telescope experiment (NITE) consists of a 16 cm rocket-borne tele- 
scope cooled by supercritical helium and a 256 x 256 low-background InSb focal plane array. 
The array is located at the focus of an //2.2 Gregorian telescope providing a 1.2° x 1.2° field 
of view with a plate scale of 17" per pixel. A 3.5 jum long- wavelength pass infrared filter 
and the intrinsic long wavelength cutoff of the InSb focal plane array defines the 3.5 to 5 
/im spectral band of the instrument. The array, sampled continuously and non-destructively 
"up the slope" at 3.8 Hz over a 10 sec reset interval, gives near-background-limited noise 
performance and maximum flexibility in data reduction. The telescope and cryogenic baf- 
fle tube were cooled by supercritical liquid helium to eliminate thermal emission from the 
instrument. The focal plane array was operated at ~ 30 K for optimal noise performance. 
Further details on the instrument may be found in Bock et al. (1998). 

NITE was designed to probe diffuse near-infrared emissions about nearby edge-on spiral 
galaxies, arising from possible low-mass stars in the massive halo. Analysis of data from two 
flights, on 1997 May 28 centered about NGC 4565, and on 1998 May 22 centered about NGC 
5907, produced strict upper limits on the near-infrared halo brightness (Uemizu et al. 1998; 
Yost et al. 2000). For purposes of assessing point source contributions to the near-infrared 
background, the two fields allow us to evaluate the relative contributions of point sources in 
the NGC 4565 (b = 86°, / = 231°; (3 = 27°, A = 177°) and the NGC 5907 (b = 51°, I = 92°; 
= 68°, A = 188°) fields. 

We observed each region by stepping the field of view 0.6° in a square pattern, and 
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integrating for 20 sec in two 10 sec exposures at each position. This sequence produces a 
complete map consisting of 8 individual exposures. The observation sequence repeats with a 
slight shift of 4' to sample the sky with different pixels. The central 32' x 36' region around 
the target edge-on galaxy is thus observed in each of the 16 exposures. Bright calibration 
stars, observed at the beginning and end of the sequence, and bright stars in the field allow 
us to calibrate the responsivity and point spread function (PSF) of the instrument. Further 
details on the observations for each flight may be found in Uemizu et al. (1998) and Yost et 
al. (2000). 



3. Data Reduction 

We present a statistical analysis from the 1997 observation about NGC 4565. The 1998 
observation of NGC 5907, at lower galactic latitude {b = 51°), suffers from a relatively larger 
contribution from stars. Furthermore, the 1998 observation has only 8 useful exposures 
due to a large terrestrial signal that appeared in the latter portion of the flight. We use 
the 1998 observations to assess the relative abundance of detected point sources at lower 
galactic latitude. The 1997 data has a noisy detector readout that samples every fourth 
column on the array. The cause of the excess noise, a loose wire-bond to the focal plane 
array, was subsequently repaired after the flight. One edge of the array was vignetted by 
a misalignment with the field stop. We carefully omit data from the noisy column and the 
vignetted region throughout the fluctuation analysis (Sections 3.2, 3.3, 3.4, and 4.2). 

We adopt the same reduction procedure as in Uemizu et al. (1998), briefly summarized 
as follows. We correct each of the 16 exposures for dark current by subtracting a master dark 
frame with a median dark current of 1.55 e _ /s, obtained by closing an internal cryogenic (5 
K) shutter just before launch. The flat field response is derived by first normalizing each 
individual exposure on the sky to its mean value. The mean is determined by applying 
an iterative 3<r clip to remove point sources. A flat field frame for each quadrant is then 
obtained by evaluating the median value in each pixel over 12 exposures. Thus 4 flat field 
frames are generated for each of the 4 pointing positions used to observe the central galaxy. 
The flat field for each pixel in this reduction procedure is thus derived from observations of 
numerous patches of sky separated by ~ 36' from the observation position. 

We divide all of the exposures by their corresponding flat field frames, and subtract a 
small constant value from each exposure to account for a slowly varying terrestrial signal. 
The time-varying signal, ~ 17 e~/s, is small compared with the total photocurrent of 160 
e~ / s. The exposures are aligned to bright stars in the field, and combined to form a 1.8° x 1.8° 
mosaic. We constructed a star mask from the mosaic, masking all pixels 3a above the average 
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over the image. Bright point sources are each masked out to regions such that their extended 
emission is no more than 2% above the total sky background based on the PSF. NGC 4565 
is similarly masked based on the results of Uemizu et al. (1998). As the combined mosaic 
has a higher sensitivity to point sources than the individual exposures, we re-evaluate the 
flat field frames using the same procedure as with the mask obtained from the mosaic. The 
new flat field frames are then combined to produce a final mosaic. 

The resulting gray-scale image is shown in Fig. 1, including data from all four array 
readouts. The regions on the corners of the image are obtained with only 2 pointings (4 
exposures), while the central region of the image is well sampled with 8 pointings (16 expo- 
sures). The image in Fig. 1 has prominent stripes due to data from the noisy column (for 
further evaluation of the statistical impact of the noisy columns see Fig 4). The remaining 
3 array readouts are separately binned and found to have very similar noise properties, with 
read noises of 52 ± 0.7 e^, 53 ± 1 e~, and 53 ± 0.7 e~ for 10 sec integration time. We 
carefully omit data from the noisy column and the vignetted regions throughout the fluc- 
tuation analysis (Sections 3.2, 3.3, 3.4, and 4.2). For the sake of simplicity, we confine the 
variation analysis in sections 3.2 and 3.3 to pixels with 100% observing time. To improve 
our statistical power, we include pixels with 75% or more observing time for the analysis in 
sections 3.4 and 4.2. 



3.1. Extraction of Point Sources 

Identification of point sources in our field of view (Fig. 1) is conducted with a PSF- 
fitting using DAOPHOT. Bright calibration stars, observed at the beginning and end of the 
flight sequence, and bright stars in the field allow us to calibrate the responsivity and PSF 
of the instrument. We iteratively extract sources by locating and masking the brightest 
sources in our field of view. The identified sources are individually masked out to regions 
where the extended PSF brightness is < 2% of the mean sky brightness, vl v = 260 nW m~ 2 
sr _1 . We use DAOPHOT to extract fluxes from identified sources. For consistency check, 
we compare the DAOPHOT fluxes with aperture photometry. In the case of very bright 
sources, we find some discrepencies, since a few are nearby extended galaxies. We thus use 
aperture photometry for the bright and extended sources. 

The resulting integrated source counts, detected with > 4a significance, are displayed in 
Fig. 2 for the 1997 and 1998 fields. The number counts are not corrected for incompleteness 
at the faint end. The large error bars at the bright end are due to the small number of 
sources. 
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3.2. Evaluation of Sky Brightness Fluctuations 

We produce two independent sky images S\ and S2 from the first 8 and second 8 expo- 
sures respectively, to discriminate sources of noise from the instrument, sky, and variations 
in the gain matrix ("flat-field" errors). We estimate the displacement between Si and S2 
accurately by comparing the centroids of bright stars. The resulting longitudinal and trans- 
verse displacements are < 0.3 and 0.2 pixel respectively. As we shall see, subtracting the 
two images allows us to distinguish noise associated with the array, compared to fluctuations 
associated with a fixed pattern on the sky. 

The measured noise in a single 10 sec integration is 53 e~, in good agreement with photon 
noise, (10s x 160e~ / s) 1 / 2 = 40e~, and readout electronics noise, ~ 30 e~ as determined 
form dark frames before launch, combined in quadrature, ((40e~) 2 + (30c - ) 2 ) 1 / 2 = 50e~. 
Instrument noise may also include various non-ideal noise sources at a low level, such as 
a fixed pattern noise on the array, low-level cosmic rays, incomplete subtraction of dark 
current, and time- varying signals. Subtracting the Si and S2 images, which are shifted in 
position relative to one another on the sky, therefore provides an accurate assessment of all 
sources of instrument noise which are necessarily correlated with array coordinates. 

To evaluate sky fluctuations statistically, we confine our analysis to the 1997 field about 
NGC 4565. Si consists of a mosaiced image from the first eight exposures, and S2 consists 
of a mosaiced image from the second eight exposures. We denote S s k y = (Si + S2) s k y /2, 
the average of Si and S2 with respect to sky coordinates, as the combined sky image. We 
mask Si and S 2 for 4a sources with a common mask derived from S s k y , taking advantage 
of the higher sensitivity to point sources in the combined image S s k y . Both Si and S 2 are 
corrected for flat-field variations as described above with a common gain matrix. We consider 
three sources of noise, instrument noise (defined to include all un-correlated noise between 
exposures), errors in the gain matrix, and fluctuations on the sky. 

Because the gain matrix is determined from observations of the sky, it can introduce 
statistical correlations in the data. We confine our analysis to the central 28' x 28' region 
to eliminate statistical correlations and to minimize gain errors, "flat-field noise" , as follows. 
The majority of the pixels in the inner region of the combined image S s k y are obtained with 
8 independent pixels on the array. The flat field for each of these pixels is obtained from 6 
different patches of sky. In constructing S s k y , each point in the inner region is compared with 
48 patches of sky located ~ 36' outside the inner region, 24 of which are independent. In the 
limit that the combined image S s k y is dominated by instrument noise, the variance of each 
pixel in the central region is increased by a factor of 8/48 = 1/6 due to flat-field variations. 
In the limit that each exposure is dominated by sky fluctuations, the variance is increased by 
a factor of 1/24. As sky fluctuations and instrument noise are actually of similar amplitude 
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in S s ky, the systematic variance of each pixel in the central region is increased by a factor 
of 1/10. The factor of 1/10 is derived by noting that flat-field errors arise from the errors 
in the off beams, with instrument variance dominating over sky variance in the off beams 
by a factor of 4. However, the variance in S sky is ~ 2 times larger than instrument noise, so 
the relative flat-field variance is (5/4)(8/48)(l/2), or approximately 1/10. Thus variations 
arising from flat field errors should generally be small in comparison to sky fluctuations in 
the inner region, and do not contribute to statistical correlations for A9 < 28'. 

We assume the variance in S s k y = (Si + S 2 ) sky /2 arises from a combination of instrument 
noise (vi), flat-field variations (vff), and sky fluctuations (v s ). For the sake of simplicity, we 
further confine the following analysis to pixels with full observing time in the central 28' x 28' 
region. The measured variance in S sky is given by 

v(S sky ) = Vi + v ff + v s = 15.04 ± 0.69 nW 2 mT A sr~ 2 . (1) 

We note that sky fluctuations will have the same amplitude in S\ and S 2 because Si 
and S 2 share a common flat field matrix. However, the gain for each pixel is derived from 
an independent arrangement of "off" positions, so flat field variations are independent from 
pixel to pixel. We assume that instrument noise is uncorrelated in the two images. The 
variance in Si and S2 is thus given as 

v(S 1 ,S 2 )=2v i + 2v ff + v 8 . (2) 

We can separate these contributions by using the fact that Si and S 2 are displaced from 
one another by 4' on the sky. By subtracting Si and S 2 with respect to sky coordinates, the 
variance in the differenced image D sky = (Si — S 2 ) sky /2 becomes 

v(D sky ) = Vi + Vff = 7.55 ± 0.35 nW 2 m~ 4 sr~ 2 . (3) 

Comparing the variance of S sky and D sky thus gives a direct measure of the sky fluctu- 
ations, v s = 7.49 ± 0.77 nW 2 m~ 4 sr~ 2 . 

We evaluate the instrument noise contribution by subtracting image pairs viewing the 
sky 10 sec apart, and determine the variance to be 185.90 ± 1.5 nW 2 m -4 sr~ 2 . As 16 
exposures comprise D sky , the estimated contribution from instrument noise is 



Vi = —(185.90 ± 1.5) nW 2 m~ A sr~ 2 = 5.81 ± 0.26 nW 2 m~ 4 sr' 2 . 



(4) 
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We then estimate the flat-field contribution to be v/f = 1.74 ± 0.44 nW 2 m~ 4 sr~ 2 . 
This result agrees within experimental error with our expectation that the variance from 
flat-field errors is Vff = (i>j + v s )/10 = 1.33, as described above. Furthermore, we can see 
that the instrument noise averages down with the number of exposures as expected. If the 
array had a fixed pattern noise that did not average down, it would be accounted as flat-field 
noise. However since the flat-field noise is in agreement with our expectation, such sources 
of excess noise must be small. The variance in each pixel of the combined image S s k y is thus 
marginally dominated by sky fluctuations, as shown in Fig. 3. 

We may test the reliability of our estimation by differencing the same two masked 
sky images with respect to array coordinates. The resulting difference image D array = 
(Si — S 2 ) array /2 is free of flat-field variation, since Si and S 2 share a common gain matrix. 
Its variance becomes 

v(D array ) = Vi + v s /2 = 10.50 ± 0.55 nW 2 m" A sr~ 2 , (5) 

in agreement with our expected variance of 9.55 ± 0.60 nW 2 m~ 4 sr~ 2 , obtained from eq. 
(1), (3), and (4) above. This consistency check shows that flat-field errors, or a non-ideal 
fixed pattern noise associated with the array, are both small and accurately estimated. It 
also shows that any effect from possible sub-pixel displacement between Si and S 2 is small 
comparing to the statistical uncertainty in our variance estimation, since eq. (5) does not 
require alignment on the sky between Si and S 2 . 

The detection of sky fluctuations can be shown graphically by evaluating the variance 
in D = (Si — S 2 )/2 as a function of the positional shift between Si and S 2 , as shown in 
Fig. 4. Pixels with 75% or more observing time were used in Fig. 4 to get more complete 
coverage, although this slightly increases i>j above the value obtained in eq. (4). The variance 
is minimized when the two images are differenced according to sky coordinates, and does 
not appreciably decrease when differenced according to array coordinates (correspondingly 
marked by the "+" in the Ax - Ay plane, Fig. 4). The width of the minimum in the variance 
is consistent with the PSF of the instrument. Therefore the observed sky fluctuations are 
unlikely to arise from improperly masking the extended emission of NGC 4565 or scattered 
light from bright stars. As a check, we changed the masking level and found no significant 
change in the above results. Also note that the striped pattern in the variance evaluation 
shown in Fig. 4, an artifact arising from omitting data obtained with noisy readout columns, 
is small (~ 5% in amplitude). 
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3.3. Auto-Correlation Analysis 

We estimate the spatial structure of the sky fluctuations using an auto-correlation analy- 
sis using the central pixels with full observing time. The two-point auto-correlation function 
(Ganga et al. 1993), 

C(0) =< 8vl v (x) ■ 5ul u (x + 9) >, (6) 

where 8vl v (x) = (ul u (x)— < vl v >), can be used to probe fluctuations on angular scales 
from 17" to 28'. We evaluate the correlation function on the summed (S sky ) and differenced 
(D sky ) images, described in section 3.1. A constant offset is removed from the summed 
image before auto-correlation analysis, obtaining an ensemble average of zero in S sky - Sky, 
instrument, and flat-field variance contribute to the variance of S s k y , whereas D s k y contains 
no sky variance, only instrument and flat-field variance. Since we only analyze the central 
28' x 28' region, we can remove any correlated structure as a result of flat-fielding. 

The resulting correlation functions are shown in Fig. 5. The correlation function at 
zero-lag, C(0), gives the variance resulting from instrument, sky fluctuation, and flat-field 
contributions as described in section 3.2. The gradual rise in Cs sky {9) near 9 = 0' is consistent 
with a combination of instrument noise and sky fluctuations arising from unresolved sources, 
spatially varying as the PSF of the instrument. Cn sky {9) shows no gradual rise, as expected 
for variance arising from instrument noise. We see no evidence for sky fluctuations at larger 
angular scales, 2' < 9 < 28'. We constrain | Cs sky {0) |< 0.5 nW 2 mr A sr~ 2 (2a), averaged 
from 10' < 9 < 28'. Our non-detection of correlated structure cannot exclude the detection 
by Matsumoto (2000) on larger angular scales (9 ~ 100') and shorter wavelengths (1.4 — 2.1 
(Km), where RMS fluctuations are significantly brighter. 

We have shown the presence of statistical fluctuations that are in excess of the instru- 
ment and flat-field noise sources. These fluctuations correlate with sky coordinates, not array 
coordinates. Furthermore the correlation analysis shows these fluctuations are spatially de- 
scribed by point sources. No known non-ideal noise source on the detector array can mimic 
these properties. The correlation analysis also shows that sky fluctuations cannot arise from 
structure in the Zodiacal foreground. 

3.4. Histogram 

Further information on the nature of the sky fluctuations may be evaluated from the 
statistical ensemble of pixels in the central region of the image. The brightness distributions 
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of the 5000 central region pixels (with 75% and more observing time) are shown in Fig. 3, 
after masking the galaxy, and detected point sources. The brightness distribution of summed 
image (Si + S2) sk /2, Hs 3ky , is noticeably skewed towards positive sources, and significantly 
wider than that of the differenced image (Si — S2) sky /2, Hu sky . The differenced histogram 
Hr) sky is well-described by a Gaussian, which eliminates noise sources on the array that could 
produce a non-Gaussian "tail" such as cosmic rays. 



4. Discussion 

4.1. Source Counts 

The source counts presented in Fig. 2 appear to be dominated by stars. The counts 
may be parameterized as 

N(S) = N (^) a , (7) 

where N = 6.0 x 10 5 sr' 1 , S = 1.0 x 10~ 6 nW m~ 2 , and a = 0.76 for the NGC 4565 
field. The index of the counts, a = 0.76, is much shallower than that of deeper .fT-band 
galaxy counts, a ~ 1.7 (Saracco et al. 1997; Bershady et al. 1998; Saracco et al. 1999; 
Jarrett et al. 2001). The extended surface brightness due to detected point sources amounts 
to vl v = 6.7 nW m~ 2 sr^ 1 in the NGC 4565 field, and 10.0 nW m~ 2 sr' 1 in the NGC 
5907 field. The ratio of the surface brightness of the two fields from detected sources, 
1.0 x 10~ 4 nW m~ 2 > vF u > 6.3 x 10~ 7 nW m~ 2 , is 1.50, somewhat larger than that 
predicted by the stellar model, 1.38, but consistent with the counts being largely stars and 
not galaxies. We are able to cross- identify 548 sources in the NGC 4565 field with Digitized 
Sky Survey images. Of the 548 cross-identified sources, 27 appear to be galaxies based on 
their measured FWHM. 

We can compare the source counts in the NGC 4565 and NGC 5907 fields to the counts 
from the model of Wright & Reese (2000) based on the Wainscoat et al. (1992) Galactic 
model, hereafter referred to as the W&R stellar model. NITE represents a significant check 
on the model, probing the counts down to vF v > 6.3 x 10~ 7 nW mT 2 (0.8 mJy) in a wide 
field. As the W&R stellar model computes differential counts in standard L and M bands, we 
note that N M (< M) = N L (< L + 0.2) over the range L = 5 to 15. Thus in order to estimate 
number counts in the NITE band, we compute N NIT e(> S) = N L (< L — 0.1), where the flux 
S is converted to nW m~ 2 in the NITE band in the same manner as the calibration stars. 
The integrated source counts agree reasonably well with the model, and are shown in Fig. 2a. 
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The integrated surface brightness, vl v {> S) — J S (dN/dS) dS, where the flux S is in vF v 
units, is compared with the measured ul u (> S) in Fig. 2b. An offset has been subtracted 
from the model by setting vl v {> S) to the measured value at vF v = 6.3 x 10~ 7 nW m~ 2 , the 
source detection limit. The offset must be subtracted to account for the surface brightness of 
bright, but rare, stars that increase the surface brightness of the model but were not present 
in the NITE field. The W&R stellar model predicts that the surface brightness contribution 
for all sources below the detection threshold amounts to Aul u = 0.7 nW mr 2 sr^ 1 for the 
NGC 4565 field. Thus we estimate that our infrared images detect 90% of the stellar surface 
brightness contribution. 



4.2. Surface Brightness Fluctuations 

The combination of a good statistical understanding of sky fluctuations and steep num- 
ber counts allows us to deeply probe the number counts of undetected sources. Following 
Condon (1974), we may recover information on the density of sources below the detection 
limit, given our large statistical ensemble of 5000 pixels (Section 3.4). We model the dis- 
tribution as arising from Gaussian instrument and flat-field noise, estimated in section 3.2, 
and unresolved sources. The surface brightness fluctuations due to Poissonian variations in 
a beam Q, where fl = 9.87 x 10~ 4 sr, 

2 1 f° ,2 dN . 

may be used to constrain the surface brightness, ul u , due to sources below the detection 
limit, 

^> s > = hf s 'w dS ' < 9 > 

Thus for a source distribution described as a power law (eq. [7]), fluctuations probe 
deeply into the background if the source counts are steep. Sources brighter than the cutoff 
limit, S c , produce a measurable change in the brightness histogram according to 

X f° C2 *V JO Q -Q.+2 

v J° Sd S 2 %dS V ' lUJ 
where Sd is the point source detection limit and S v is the change in the variance produced by 
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all sources fainter than S c . Since our measurement of the fractional uncertainty in the sky 
variance is 5 v /v = 10.3% ( v s = (7.49 ± 0.77) nW 2 m~ A sr~ 2 , Section 3.2), the histogram is 
sensitive to sources much fainter than the detection limit Sd- For example S c /Sd = 10~ 2 in 
the case of Euclidean counts (a = 1.5). If a > 2, fluctuations will be dominated by sources 
where the power law relationship eventually rolls over, in principle well below Sd- 

The width and near-Gaussian shape of the Hs sky in Fig. 3 indicate a rapid increase in 
the source counts below the detection threshold. Comparison with i^-band galaxy counts 
(a ~ 1.7) also suggests the need for higher a. Therefore, we fit the data to two components, 
stellar counts taken from the W&R stellar model, and a power law population (eq. [7]) 
with undetermined N and a. We assume that the latter component represents galaxies. 
Assuming the source counts are dominated by stars taken from the W&R stellar model, our 
simulation results in a poor approximation to the observed fluctuations, as shown by the 
solid line in Fig. 6. 

We carry out a P(D) analysis to estimate the density of sources below the source detec- 
tion limit, and constrain the associated background produced by these sources. We created 
a Monte Carlo sky to best estimate faint galaxy source count parameters N and a (eq. 
[7]). The simulated sky contains point sources with flux 1.0 x 10~ 9 nW m~ 2 < vF v < 
6.3 x 10 -7 nW m~ 2 . According to eq. (10), we must include all sources brighter than S c . 
Including sources fainter than S c has no appreciable affect on the results. We begin our 
simulation with a 1° x 1° field, gridded with 4.25" pixels. We populate the grid with ran- 
domly placed delta functions, which we generate assuming a galaxy power law source counts 
(eq. [7]), and the W&R stellar model, selected in steps of 0.1 magnitude (corresponding to 
a multiplicative factor of 1CT - 04 nW mr 2 ). We then convolve the simulated image with the 
PSF, and re-bin to 17" pixels to match that of the detector array. Finally we add Gaussian 
noise with dispersion a n to each pixel to simulate instrument and flat-field noise. 

The simulated sky is then compared with the observed histograms Hs sky and Ho sky (see 
Fig. 3). We fit the histograms jointly, allowing N , a, 1$, and a n to be free parameters. 
I accounts for the difference in mean between the simulated and the observed skies. We 
must let I be undetermined since our data contains a significant undetermined contribution 
from Zodiacal emission. Letting I vary also means that our simulation is insensitive to the 
faintest sources (S < S c ) applied in the Monte-Carlo simulation, since these sources add 
surface brightness but give no appreciable change in the histogram. We let a n vary to allow 
for uncertainty in the instrument and flat-field noise. The quantity a n has an experimental 
error (see eq. [3]) that affects both H Ssky and H Dsky . Fitting Hs sky and H Dsky jointly properly 
accounts for the uncertainty in a n , and essentially constrains a n to values within experimental 
errors. 
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We determine the best fit (reduced x 2 ) to the histograms for a given set of N and a 
values, marginalizing over a n and I . We semi-analytically truncate the analysis to a < 1.8, 
based on eq. (10), to handle the case where the histogram is sensitive to sources fainter than 
our minimum Monte-Carlo source flux. We repeat simulations multiple times to accumulate 
statistical power and to reduce the numerical error in the simulation. We examine the 
goodness-of-fit over a wide range of N and a, as shown in Fig. 7, with 50 simulations 
performed for each point in the N — a parameter space. Our best-fit result, as shown by 
the dashed line in Fig. 6, and marked by the "+" in Fig. 7, gives N = 0.8^6 x 10 5 sr^ 1 , 
and a = 1. 65^35, with So fixed to the numerically convenient value of 10~ 6 nW m~ 2 . 

The error ellipse shown in Fig. 7 may be translated into iV(> S) and vl v (> S) vs. S 
plots by applying eq. (7) and (10), and propagating the errors appropriately. We evaluate 
the "true" cutoff fluxes using a Monte-Carlo simulation, and find them to be brighter than 
expected from eq. (10), presumably due to the addition of the fitted parameters Io and a n . 
We adopt the simulated cutoff fluxes for the la and 2a error contours, rather than those 
from eq. (10). We truncate the error bars when any portion of the error ellipse falls below 
the corresponding cutoff flux. The resulting constraints are shown in Fig. 8, along with the 
source count data from the NGC 4565 field. The associated errors on N(> S) and vl v {> S), 
also shown in Fig. 8, improve at fainter fluxes because the error ellipse rotates to make the 
projected error in N(> S) and vl„(> S) smaller. We thus determine the integrated source 
counts N(> S) = 1.79+^ x 10 7 sr~\ contributing vl v {> S) = lM±° f 4 nW m~ 2 sr' 1 to 
the CIBR, evaluated at S = 4.0 x 10" 8 nW m' 2 . 

We note that this Monte-Carlo technique is limited by the statistical power of the data. 
A more precise statistical determination of the histograms Hs 3ky and Hi> sky would produce 
a smaller error ellipse in Fig. 7. This would allow us to exclude the smallest values of a 
currently in the error ellipse, and consequently carry the curves in Fig. 8 deeper since the 
maximum allowed cutoff flux would be lower. 

The values of iVo and a agree with the observed galaxy counts at fT-band (N(> Sk) ~ 
1.8 x 10 7 sr- 1 at S K ~ 1.0 x 10~ 7 nW m~ 2 , and a ~ 1.7; Saracco et al. 1997; Bershady et 
al. 1998; Saracco et al. 1999; Jarrett et al. 2001), giving an average color K — M NITE = 1.0. 
The i^-band counts show a roll-off at K ~ 17, which is not included in our source count 
model. Alternately, we can assume reddened if-band galaxy counts (K — M^ite = 1-0) in 
our simulation. The result is consistent with the histogram data and produces a background 
of 2 nW mT 2 sr^ 1 . 

As a check, we fit to the histogram Hg and employ three point Gauss-Hermite quadra- 
ture summation (Zwillinger 1996; Ganga et al. 1997) to estimate the errors arising from 
variation in a n . The best-fit result is consistent with the joint analysis. 
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4.3. Limits on Extragalactic Background 

DIRBE produced an absolutely calibrated all sky map at 10 different IR wavelength 
from 1.25 fim to 240 fim. Unfortunately, DIRBE's 0.7° beam is contaminated by Galactic 
foreground stars. Since it is difficult to observe large regions of sky to sufficient depth 
in the thermal near-infrared, models and fits of the Galactic stellar flux have been the 
primary method used to remove the contamination from the DIRBE maps (Hauser et al. 
1998; Dwek and Arendt 1998; Wright and Reese 2000). Recently the Two Micron All Sky 
Survey (2MASS) provided the needed large-scale, high-resolution data to help determine the 
background at 1.25 fim and 2.2 fim (Cambresy et al. 2001; Gorjian et al. 2000). Gorjian, 
Wright, & Chary (2000) also determined the brightness contribution from foreground stars 
at A = 3.5 fim in selected DIRBE fields, down to a limiting flux of 3.8 x 10~ 5 nW m~ 2 . In 
comparison, these observations are 50 times deeper in vF v than Gorjian, Wright, & Chary 
(2000) and cover a wider field. 

Galactic dust emission appears to be negligible in our field. Arendt et al. (1998) estab- 
lished an empirical relation between residual Galactic dust 100 fim emission not associated 
with neutral hydrogen, and emission at other wavelengths. We use 16 DIRBE pixels on 
the perimeter of our field to derive a 100 fim intensity of 52.12 nW mr 2 sr" 1 , of which 
32.32 nW mr 2 sr~ l remains after extrapolating to N HI = 0. The inferred contribution 
from dust emission is completely negligible, 0.06 nW mr 2 sr^ 1 at 3.5 fim and 0.09 nW mr 2 
sr^ 1 at 4.9 fim. Likewise the emission from the edge-on spiral galaxy in the center of our 
field makes a small contribution to the surface brightness. Although technically part of the 
extragalactic background, we exclude the flux, vF v = 7.4 x 10~ 4 nW m~ 2 for NGC 4565 
and vF v = 4.2 x 10~ 4 nW m~ 2 for NGC 5907, amounting to 0.75 and 0.43 nW m~ 2 sr' 1 
averaged over the 1.8° x 1.8° field, respectively, from the CIRB. 

The image obtained by NITE provides us with a large scale, high-resolution image, 
and so we can use it to subtract stellar flux from a DIRBE map with zodiacal subtraction 
(Kelsall et al. 1998). We align the NITE and DIRBE fields, and use the nine DIRBE pixels 
completely contained within the NITE image to calculate the CIRB. 

In Table 1 we account for the stellar foreground in our two fields. We remove the small 
flux from galaxies resolved in POSS plates from the stellar foreground. The color of the 
stellar foreground produced by 5 < L < 13 stars is L — M = —0.12 in the Wright & Reese 
(2000) model (note the value for vl v {> S) is slightly different than for iV(> S)). Undetected 
stars, as noted in sections 4.1 and 4.2, increase the stellar surface brightness by only 10% 
based on the stellar model. Assuming this color, the total estimated stellar foreground is 
11.6 ± 0.5 (4.1 ± 0.2) nW mr 2 sr' 1 at 3.5 (4.9) ^m for the NGC 4565 field, and 17.4 ± 0.8 
(6.2 ± 0.3) nW mr 2 sr' 1 at 3.5 (4.9) fim for the NGC 5907 field. 
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In Table 1 we estimate the residual emission, subtracting the stellar foreground from 
the zodi-subtracted DIRBE sky brightness in our fields. We also account for the small 
contributions from detected galaxies, undetected stars, and the central edge-on galaxy. The 
resulting residuals are 24.0 ± 2.7 (36.8 ± 6.0) nW m~ 2 sr' 1 at 3.5 (4.9) ^m for the NGC 
4565 field, and 14.0 ± 2.6 (22.2 ± 5.9) nW mr 2 sr' 1 at 3.5 (4.9) ^m for the NGC 5907 
field. The errors are dominated by the quoted zero-point errors from the zodiacal subtraction 
(Kelsall et al. 1998). 

The large discrepancy between the NGC 4565 and NGC 5907 fields may arise from 
the zodiacal subtraction. The NGC 4565 field lies at low ecliptic latitude (j3 = 27°) where 
zodiacal residuals are more pronounced. We confirm that a bright star, prominent in the 
DIRBE map and just outside of the NGC 4565 field, does not contribute significant emission 
to the DIRBE pixels we used to estimate the sky brightness. 

The NGC 5907 residual agrees with residuals reported by several authors: Wright & 
Reese (2000), vl v = 12 ± 3 nW mr 2 sr~ l at 3.5 ^m (CIRB detection), Gorjian, Wright & 
Chary (2000), vl v = 11.0 ± 3.3 nW rrT 2 sr' 1 at 3.5 /im (tentative CIRB detection), and 
Dwek & Arendt (1998), vl v = 14±3 nW m' 2 sr' 1 at 3.5 /im, and vl v = 25±6 nW mr 2 sr' 1 
at 4.9 fim (tentative CIRB detection and upper limit, respectively, with updated values taken 
from Hauser & Dwek 2001). The NGC 5907 residual also agrees with the residual brightness 
in the DIRBE HQB regions, vl v — 11 ± 6 nW m~ 2 sr^ 1 at 3.5 /im and vl v = 25 ± 8 nW 
m -2 gr -i a -j. ^ g k u ^. jj auser e ^ a j 1998 report only upper limits on the CIRB. Although 
the NGC 5907 residual adds to the circumstantial evidence for a detection of the CIRB, we 
cannot test for isotropy in our small field. 

We estimate the CIRB by accounting for the sources producing the observed surface 
brightness fluctuations. Because we assume a truncated power law for the source counts, we 
can only constrain sources down to S c , and thus produce a lower limit on the CIRB. The 
rapidly rising source counts required to explain the histogram are very clearly due to galaxies, 
as expected from K-band survey counts. We evaluate the detectable infrared background 
by integrating over sources down to the simulated la cutoff flux, S c = 4.0 x 10~ 8 nW mT 2 , 
16.7 times below the point source detection flux, contributing vl v {> S) = 1. 041^34 nW mr 2 
sr" 1 to the cosmic infrared background. 

5. Conclusions 

We analyze the point source contribution to deep, wide-field 4 jum images. Detected 
sources are dominated by stars, and are in general agreement with predictions from the 
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Wright & Reese (2000) stellar model. We report a detection of surface brightness fluctuations 
detected with an ensemble of 5000 pixels. The brightness fluctuations do not show correlated 
structure and appear to be dominated by galaxies. We estimate the C1RB arising from 
sources that contribute to the fluctuations. SIRTF can use sky fluctuations to study the 
cosmic infrared background in greater detail, probing galaxy counts to fainter limits. Shallow 
SIRTF images, distributed over the sky, can remove the stellar foreground from D1RBE, 
allowing a test of the isotropy of the residuals. 
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Fig. 1. — Gray scale image of the full 1.8° x 1.8° field centered about the edge-on spiral 
galaxy NGC 4565. The image is rotated so that equatorial north is up and east to the left. 
Data from the noisy readout and the vignetted region are omitted in the fluctuation analysis 
(Sections 3.2, 3.3, 3.4, and 4.2). 



- 21 - 




Fig. 2. — a) Integrated number counts are plotted for all sources detected with > 4<r sig- 
nificance, shown for the NGC 4565 (a) and NGC 5907 (o) fields. The solid curves denote 
predictions of the Wright and Reese (2000) stellar model, b) Integrated surface brightness 
is plotted for the two fields, with data and model predictions as in Fig. 2a. 
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Fig. 3. — Brightness distributions of masked IR image about NGC 4565. The histogram plot 
of the summed image (Si + S2) sky /2, Hs aky , exhibits an asymmetry about the mean, and is 
significantly wider than that of the differenced image (Si — S 2 ) sky /2, H Dsky . The difference 
histogram H Dgky is fitted to a Gaussian. The variance in each pixel of the masked image is 
thus marginally dominated by sky fluctuations. The near-Gaussian shape of H Ssky indicates 
a large number of sources per pixel, requiring a rapid increase in the source counts below 
the point source detection threshold. 
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Fig. 4. — Variance in D = (S± — S-i)/2 is plotted as a function of positional displacement 
between S\ and S2. The Ax - Ay coordinates correspond to the relative offsets on the sky 
of Si and S 2 ■ The minimum at (Ax, Ay) = (0, 0) pixel corresponds to subtracting the 
images in sky coordinates, confirming the presence of a fixed pattern on the sky, or "sky 
fluctuations". The shape of the minimum is consistent with the instrument PSF FWHM. 
There is no appreciable decrease in variance at (Ax, Ay) = (—12, —1), when the two images 
are differenced according to array coordinates, as marked by the "+" in the Ax - Ay plane. 
The overall fluctuation (~ 5% in amplitude) pattern in variance evaluation is an artifact 
arising from omitting data obtained with noisy readout columns. 
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Fig. 5. — Auto-correlation function of the sum and difference frames. The gradual rise in 
the sum frame, Cs sky {9) near 9 = 0', is consistent with a combination of instrument noise 
and sky fluctuations arising from unresolved sources, spatially varying as the instrument 
PSF. Co sky (9) shows no gradual rise, as expected for variance arising from instrument and 
flat-field noise. 
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Fig. 6. — Data points with Poissonian error bars denote the measured histogram Hs aky of 
the combined sky image. The estimated histogram for the combination of instrument noise 
and undetected stars taken from the W&R stellar model results in a poor approximation to 
the observed fluctuations, as shown by the solid line. Additional faint sources are required to 
broaden the simulated histogram. We model the observed fluctuations with two components, 
stellar counts taken from the W&R model, and a power law galaxy population as described 
in section 4.2. The best-fit simulated distribution is shown by the dashed line. 



-26- 




a 



Fig. 7. — 68% and 95% confindence contours for our best-fit result to the source population 
producing brightness fluctuations, as described in section 4.2. Our best-fit result, marked 
by "+", gvies iVo = x 10 5 sr" 1 , and a = 1.65t{j:S|j, with S held at the numerically 

convenient value of 10 -6 nW m~ 2 . 
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Fig. 8. — Estimated N(> S) and vl u (> S) curves based on the Monte- Carlo P{D) analysis 
derived from the error ellipses shown in Fig. 7. The solid curves denote the best fit result, 
with dashed and dotted curves indicating 68% and 90% confidence curves respectively. The 
curves are terminated at S c determined from the Monte-Carlo simulation. Since the la and 
lcr error ellipses include low values of alpha, with a higher corresponding cutoff flux, these 
curves are truncated at their corresponding maximum cutoff flux which is higher than the 
cutoff flux for the best fit curve. Also plotted are the data points from sources extracted 
from the NGC 4565 field, the field used to produce the histograms for the P{D) analysis. 
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Table 1. Estimation of the Sky Brightness and Foreground Components 



Component 


3.5 (im 


NGC 4565 

3.5-5 


4.9 /j,m 


3.5 (im 


NGC 5907 

3.5-5 /im 


4.9 (im 


DIRBE (Zodi sub.) 


37.0 1 ± (1.7)(2.1) 2 




41.5 ± (0.8)(5.9) 


32.1 ± (1.3)(2.1) 




28.7 ± (0.5) 


ISM 


0.06 




0.09 


0.04 




0.06 


Detected Stars 


10.5 ± 0.5 


6.7 ± 0.3 


3.7 ± 0.2 


15.7 ± 0.8 


10.0 ± 0.5 


5.6 ± 0.3 


Undetected Stars 


1.1 




0.4 


1.7 




0.6 


Detected Galaxies 


0.20 


0.13 


0.07 








Galaxy 3 


1.18 


0.75 


0.41 


0.66 


0.42 


0.23 


R csidual 


24.0 ± 2.7 




36.8 ± 6.0 


14.0 ± 2.6 




22.2 ± 5.9 



1 The units for all values are in nW m 2 sr 1 

2 The first number in the parentheses is the random la error and the second one is the systematic la error due to zodiacal 
subtraction from (Kclsall et al. 1998) 

3 The NITE flux was converted to the L and M band fluxes assuming the same colors as those for the undetected stars in the 
Milky Way. 



